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^ i Abstract. 

^! We are concerned with large scale magnetic field dynamo generation and propagation 

O 
U 



of magnetic fronts in turbulent electrically conducting fluids. An effective equation for the 
large scale magnetic field is developed here that takes into account the finite correlation 



^ . times of the turbulent flow. This equation involves the memory integrals corresponding 

to the dynamo source term describing the alpha-effect and turbulent transport of mag- 
netic field. We find that the memory effects can drastically change the dynamo growth 
rate, in particular, non-local turbulent transport might increase the growth rate several 
times compared to the conventional gradient transport expression. Moreover, the inte- 
gral turbulent transport term leads to a large decrease of the speed of magnetic front 
propagation. 
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The problem of the generation and propagation of a magnetic field in turbulent electri- 
cally conducting fluids is of fundamental importance due to various applications in plasma 
physics, astrophysics, geophysics, etc. It has attracted much attention since the 

late 60s when it was realized that the conditions for the occurrence of a large scale mag- 
netic field can be found by applying a simple technique involving the mean helicity of the 
turbulence [^f. This technique is based on the effective macroscopic equation governing 
the large scale magnetic field B. The standard form of the mean field dynamo equation 
is 



O'B 

— = rot(aB) + f3AB + rot[u x B] (1) 

where u is the mean velocity field, a is the helicity and (3 is the turbulent diffusivity. 
This equation is a common starting point for analyzing the generation of the large-scale 
magnetic field pl- |H|- It has been also used for the analysis of the propagation of magnetic 



fronts in spiral galaxies H- |[Tl 



The main disadvantage of equation (1) is that it has been derived for turbulent flow 
involving only two separated length scales for the velocity field - the integral length scale 
and the small turbulent scale . It is clear that the assumption of two separated scales is 
rather unrealistic for fully developed turbulent flow, which involves a continuous range of 
spatial and temporal scales [^. Thus, the purely local equation (1) is applicable only under 
the assumption of a clear-cut separation between macroscopic behavior of the averaged 
magnetic field and the turbulent fluctuations at the "microscopic level". 

It should be noted that the equation (1) is similar to the convection-diffusion-reaction 
equations |T^. In fact, it can be reduced to a famous Fisher-Kolmogorov-Petrovskii- 
Piskunov (FKPP) equation |^ which has become a basic mathematical tool in the 
theory of propagating fronts travelling into the unstable state of the reaction-diffusion 
systems. There has been an increased interest in this topic, because of the large number 



of physical, chemical and biological problems that can be treated in terms of the FKPP 
equation (see, for example, [jl2[-|l^). Recently, there has been a tremendous activity to 



extend this analysis by introducing more realistic description of the transport processes. 
The main motivation for this is that the diffusion approximation for transport admits an 
infinite speed of propagation. Due to this nonphysical property of the diffusion solution, 
the FKPP equation yields an overestimation of the propagation speed of travelling fronts 
14|-[1^. We expect that a similar situation might take place in magnetohydrodynamics. 



The mean-field dynamo equation (1) also admits an infinite speed of a magnetic field 
propagation. It is clearly a nonphysical property, because the speed of magnetic field 
propagation can not exceed the velocity of the largest eddy of turbulent flow. The origin 
of this contradiction lies in the 5— correlated-in-time approximations for the turbulent 
velocity field (see below). In reality these correlations have finite times of relaxation, and 
what is more these might be of the same order as the characteristic times for the growth 
rate of the large scale magnetic field since the physical origin of these correlations and the 
magnetic field generation are the same, namely, the turbulent fluctuations. 

It is the aim of this Letter to extend mean field dynamo equation (1) to the case when 
long range in time correlations of turbulent flow are taken into account and find out how 
the non-local in time effect might influence the critical conditions for the generation of 
magnetic field and its spatial propagation. 

The most general phenomenological formulation of the dynamo problem that is con- 
sidered here is represented by the equations for the average magnetic field B 0: 

— = loiS + rot[u X B] , = (u' X B') , (2) 

where primes denote the turbulent fluctuations and the angular brackets denote an ensem- 
ble averaging over the turbulent pulsations. The main closure problem here is to express 
the turbulent electromotive force £ in terms of the average field B. The classical expres- 
sion leading to (1) is £^ = a B — /? rotB, where a and (3 are the statistical characteristics 



of the turbulent velocity u'. However, under the assumptions of infinite conductivity and 
weak turbulence the electromotive force £ can be written as IHl: 



1 /•* 

£^(x, t) = — / (u'(x, t)- rot u'(x, s)) B(x, s) (is — 

3 J -oo 

1 /•* 

— / (u'(x, t) ■ u'(x, s)) rotB(x, s) ds . 

3 J-OO 



(3) 



The local mean field equation (1) can be derived from (3) under the assumptions 
that the correlations appearing in (3) are approximated by the delta-functions in time: 
(u'(x, t) ■ rot u'(x, s) = —3a 6{t — s), (u'(x, t) ■ u'(x, s)) = 3f3 6{t — s). However, as was 
mentioned in [Q (p. 136) "... the assumption of instantaneous correlations seems to be 
a serious restriction to the theory...", since in real turbulence the characteristic times of 
these correlations are finite. It follows from (3) that in the case of finite correlation times 
the electromotive force should contain the integrals over the history of B and rotB. 

In this Letter we suggest the following general form for the electromotive force S in 
the limit of infinite conductivity 



8 = a(x) J ( V^) ~ ^^^^ / ^'^ ( V^) ^) ' 

where Ga iy) and Gp (y) are positive, decreasing functions that tend to zero as ?/ — > oo. 
The parameters and control the time correlations in the random velocity field at 
a fixed space position. By inserting (4) into (3) one can get the non-local in time mean 
field dynamo equation: 



^ = f_ rot [aBis)] ds - (5) 

j Gf3 (——^ [/^ rotB(s)] ds + rot[u x B] . 
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It should be noted that this equation is vahd only in the limit of infinitely large conduc- 
tivity. This case is of specific interest for magnetohydrodynamics of plasma and many 
problems of astrophysics and geophysics [^. Many different constitutive models might 
arise from different choices of (y) and G/3 (y) . The equation (1) can be considered as a 
limiting case of (5) when Ga^i3{t — s) = 6{t — s) (tq^^ 0) . The transport memory kernel 
Gi3 ensures the finite velocity of propagation of a magnetic field, which is determined by 
the rate of turbulent pulsations 0. The relaxation time rg can be determined from the 
root mean square velocity Uj-ms = \fWF^p of turbulent fiuctuations. The integral kernel 
GoL is used to express the fact that the dynamo growth of a magnetic field at a time t 
in the local vicinity of the spatial point x is determined by the past values of rotB. 

Although there has been a large number of studies of the generation of macroscopic 
magnetic fields, the implications of finite correlations for the generation and propagation 
of the large scale magnetic field have not been studied. The implications of our results for 
dynamo theory are twofold. First, the non- locality in time greatly influences the dynamo 
growth rate of a magnetic field, however, the generation conditions are not sensitive to 
the addition of the memory effect due to non-zero time correlations: they are robust to 
the addition of non-local terms for the mean field dynamo equation (5). Secondly, the 
memory effects introduce drastic changes to the front dynamics of a magnetic field. 

In what follows, we study the influence of the memory effects on the dynamo generation 
and propagation by using an important example of a thin turbulent slab of thickness 2/i 
and radius R (i? ^ /;,), which rotates with the angular velocity loir). This is a standard 
model for disc-like galaxies. We neglect the effects of compressibility, diamagnetism and 
deviations from the axial symmetry. We use the standard approximation of constant 
turbulent diffusion coefficient /? ^, |], [ill . We restrict our analysis to the kinematical 



aspects of the problem, neglecting the influence of the magnetic fleld on the turbulent 
flow, i.e. the dependencies of a,/?, Tq and on the magnetic fleld B. With these 
simpliflcations, the equations for the mean axisymmetric magnetic fleld in polar cylindrical 



coordinates (r, y9, z) with z-axis coincident with the rotation axis follow from the basic 
equation (5): 



dBr 
~dt 



(6) 



G, 



H 

Qj. 



r or 



ds 



(7) 



Gj3 
oo V ''73 



t-s\ (d'^B,„ d 



+ 



Q^2 Qj, 



r or 



ds 



Here g = rduo/dr is the measure of differential rotation, and we are interested only 
in Br and B^p components of a magnetic field, because Bz/B^^ip = 0{h/R) 0, |^. 
These components obey the vacuum boundary conditions on the thin disc surfaces [|, |^: 
Br,Lp{z = ±h) = 0. The main goal of our analysis is to find the local growth rate of the 
magnetic field and its spatial propagation, taking into account the memory effects. Thus, 
considering the field generation in the usual way ||, |], |10|, we neglect first the radial 
derivatives in equations (6), (7). We represent the components of the magnetic field as 
follows Br,yy{z,t) = br,^{z) expi^j t), where br,<p{z) have to be found from the eigenvalue 
problem 



(7 +^)br = -Rc.^^ , ((f+^)K = R^br + Ra^^ , = ±1) = , 

(8) 
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where 7 = 7/ fp{lTp) , = R^/ fpi'^Tp) , = Rafa{lTa/ Ra)/fp{lTp) , faAl^afi) = 
/o°° G'q:,/3(^/^q:,/3) exp(— 7^) , Oa = T^/Ra , = ■ Here we use the dimension- 
less variables z — > z/h,t (3t/h?^a —>■ aoa{z), and dimensional parameters = 
aQTa/h,Tj3 = (3Tp/h?,Ra = a^h/ [3,Ryj = gh'^ / (3. The parameter 7 describes the growth 
rate of the magnetic field. The eigenvalue problem (8) reduces to the well-known form 
0, 1^, |lOl in the case when T^,/? = 0. Otherwise, the eigenvalue 7 becomes a function 



of these correlation times, and this function is determined from the problem (8) using the 
renormalized parameters ^,Ra,Rw 

The problem (8) coincides with the generation equations in the local mean-field dy- 
namo theory H, ^ |l^. The only difference is that the renormalized parameters j,Ra 
and Ry^ are dependent on the increment 7 of the time growth of the magnetic field. This 
means that the generation equations (8) are universal and do not depend on the specific 
choice of the memory kernels Ga,i3- The main physical conclusion from this universality is 
that since physically meaningful memory kernels must satisfy the normalization condition 
fa^pil = 0) = 1, then the threshold combinations of dynamo parameters Ra and Rp 
at a given a{z), providing the critical point of instability (Re7 = 0), are the same for all 
memory kernels. These critical parameters can be determined from the local mean-field 
dynamo theory. To determine the growth rate 7 in the generation region, the eigenvalue 
problem (8) should be solved as 7 = ^{Ra^Rw)- This relation with the definitions of 
the renormalized parameters 'j,Ra,Rw leads to a transcendental equation for 7. The 
specific dependence of 7 on the relaxation times Tq,^^ and dynamo parameters Ra,w 
has to be determined by the specific forms of the memory kernels Ga,i3- 

Let us illustrate the general results by using the important example of auj dynamo 
{Ra <^ l-Rtol), which is of specific interest for astrophysical problems f^l-iUj We 
use the exponential relaxation form 

GaAv) = exp {-y I Ta,p) I Ta,p . (9) 
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The solution of the eigenvalue problem (8) depends on the form of the function a{z). The 
au approximation requires the anti-symmetric character a{z) = —a{—z). Following the 
asymptotic method of solution of the eigenvalue problem (8) ^ [1^, for the model case 
a{z) = z, the increment of maximal growth satisfies the algebraic equation: 



The growth rate 7 is shown in Fig. 1 as a function of dimensionless correlation times 
Tq, and Tj3. Clearly, the non-local terms corresponding to the turbulent helicity and 
the turbulent transport influence the growth rate 7 in opposing way. The turbulent a- 
source non-locality leads to a slower field growth (curve 1) due to the physical nature of the 
process: the integral representation of the a-term in basic equations (4), (5) introduces the 
time memory of excitation. The non-locality of turbulent /5-transport reduces the energy 
loss out of the disc and, hence, results in an increase of the field growth rate (curve 2). 
Moreover, for dynamo numbers D, close to the critical value Dcr, the influence of the 
transport non- locality becomes very significant, the growth rate 7 increases drastically 
with Tf^. It is worth noting that the critical value Dcr — of the dynamo number 

D, at which the generation commences (7 = 0), does not depend on the memory effects; 
this result is in agreement with the general conclusion made above. 

Now we are in a position to discuss the problem of magnetic front propagation. If the 
dynamo excitation takes place within a certain radius (say, r < tq) then the magnetic 
field can propagate in the form of travelling fronts According to the local dynamo 
theory the speed f of a propagating magnetic front is proportional to y/J^y due to the 
FKPP type of the mean field equation (1). The problem is that in the local mean-field 
description of magnetic fields, turbulent diffusion induces an infinitely long Gaussian tail 
ahead of magnetic front which leads to the overestimation of the propagation rate. A 




D 



■w 



(10) 
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qualitatively different situation appears when the memory effects are fallen into account. 
This means that the speed v can not exceed the maximum possible velocity \J (^/Tp 
however large the growth rate 7 is. 

Now let us turn to the problem of front propagation in the case when the dynamo 
equation with memory (5) is considered in the thin disk approximation. Assume that the 
initial distribution of magnetic field satisfies: B = Bq = const if r < tq and B = 
otherwise. The main quantity of interest is the speed v at which magnetic field propagates 
in the form of a self-similar wave for the large values of r and t . In the large-distance 
limit (r 00) the radial Laplace operator ^(^^'") in equations (6) and (7) can be 
approximated by the second derivative ^ and we may find the solution of these equations 
in the form of Fourier modes Br,^{r, z, t) = K^^^z) exp(7 t) exp{ikr). Substitution of this 
expression into equations (6) and (7) leads to the problem for eigenfunctions br^^ IQ?!!]) 
from which one can find the equation for the exponent '~f{k) as a function of the wave 
number k. The general theory of front propagation in non-local reaction-diffusion media 



Iq , |16|, based on the saddle-point method of calculation of the inverse Fourier integrals, 
leads to the following equation for the propagating velocity v = 7(A)/A , where X = ik 
satisfies the equation 7(A)/A = dj{X)/dX. We found that the front velocity v is 
a monotonically decreasing function of both correlation times Tq^^ . In particular, it 
depends weakly on Tq, (see Fig. 2) , but decreases significantly (up to 2-3 times) with 
growing . Therefore, the use of FKPP-like estimation for the travelling wave velocity 
V ~ VJ3^ for systems with memory leads to significant overestimations. This general 
conclusion is of great importance not only for magnetic dynamo front propagation but 
for a wide class of excitable media. 

Basically, we have extended the classical mean field dynamo theory to the case when 
the memory effects are taken into account. We have suggested the effective equation for 
the large scale magnetic field that takes into account the finite correlation times of the 
turbulent flow. This equation involves memory integrals corresponding to the dynamo 

9 



source term describing the alpha-effect and turbulent transport of magnetic field. We have 
found that memory effects can drastically change the dynamo growth rate, in particular, 
the finite turbulent transport involving memory might increase the growth rate several 
times. We have also found that memory effects lead to the essential decrease of the speed 
of magnetic front propagation. 

The authors thank David Moss for interesting discussions. This work was carried out 
under EPSRC Grant GR/M72241 and with the financial support of the CRDF Award 
No. REC-005. One of the authors (AI) was partly supported by the RME Grant No. 
EOO-3.2-210. 
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Figure 1: Dependencies of the dimensionless growth rate 7 on the turbulent a-source 
relaxation time at T/3 = (curve 1) as well as on the transport time T/j when Tq, = 
(curve 2). The memory kernels Ga^p are chosen to have form (9); the dimensionless 
dynamo parameters of the system are: Ra = 1, Rw = 6.1. 
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Figure 2: Dependence of the dimensionless front velocity v = v(Ta,Tj3)/v{0,0) on the 
relaxation time T/3 at = (curve 1); 0.3 (curve 2) and 0.9 (curve 3). The dynamo 
parameters are Ra = 2 and R^ = —40 
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